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Abstract 

In this paper, we derived a critical condition for matter equilibration in heavy ion 

collisions using a holographic approach. A gravitational shock waves with infinite trans- 
it , 

verse extension is used to model infinite nucleus. We constructed the trapped surface 
in the collision of two asymmetric planar shock waves with sources at different depth 
in the bulk AdS and formulated a critical condition for matter equilibration in collision 
of "nucleus" in the dual gauge theory. We found the critical condition is insensitive to 
the depth of the source closer to the AdS boundary. To understand the origin of the 
critical condition, we computed the Next to Leading Order stress tensor in the bound- 
ary field theory due to the interaction of the nucleus and found the critical condition 
corresponds to the breaking down of the perturbative expansion. We indeed expect 
non-perturbative effects be needed to describe black hole formation. 
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1 Introduction 



The AdS/CFT correspondence is conjectured as a duality between weakly coupled grav- 
ity theory and strongly coupled AA=4 Super Yang-Mills theory in the limit of large iV c 
and strong couplingfl, 3, 2]. Its applications to strongly coupled Quark Gluon Plasma 
(sQGP) have revealed many novel features of the strongly coupled medium such as very 
low viscosity [4], absence of jets, Mach cone formation and other hydrodynamical phenom- 
ena, see e.g. review [5]. While static and near equilibrium properties have been extensively 
studied in this context, deriving corrections to hydrodynamics, the out-of-equilibrium as- 
pects of the strongly coupled gauge theory remains less understood. One of the main 
challenges in heavy ion collisions remains the understanding of early equilibration of matter 
produced in the collisions. 

Recently there have been several attempts to model the initial non-equilibrium stage 
of the collision, including [6] and [7]. Our paper [8] provided relatively simple description 
of black hole formation, due to elastic membrane falling under its own weight. In this case 
the ultimate equilibration is always assured, and it happens as a gradual propagation of the 
equilibration boundary in the scale space (along the 5-th holographic dimension), from the 
ultraviolet (UV) toward the infrared (IR) direction. 

According to principles of AdS/CFT, due to large N c limit all issues have to be 
understood in terms of classical gravity problem. Thermal equilibration of matter and early 
entropy production is in this setting dual to the formation of a (black hole) horizon, trapping 
some amount of information from the distant observer, where our world is. This mechanism 
not only is able to provide some lower bound on the amount of entropy production in the 
collision, but it also provides qualitative "yes" or "no" answer if the information trapping 
does or does not happen, as a function of given initial condition of the problem. Thus 
one of the interesting unexpected features of the problem are some rapid transition into a 
new regime, as a function of e.g. collision energy, density of the colliding objects or (not 
discussed in this work) the impact parameter of the collisions. 

The relation between the trapped surface at the collision moment and the lower 
bound on the entropy production has been introduced by Gubser, Pufu and Yarom [9], 
who have considered collision of ultrarelativistic small black holes in AdSs. It can be 
viewed as a collision of gravitational shock waves, having near-zero longitudinal width but 
possessing a certain profile in 3 transverse coordinates x 2 , x 3 and the holographic coordinate 
z. Mathematically, the trapped surface at the collision point satisfies the Laplace eqn, 
plus certain nontrivial boundary conditions on the surface. Furthermore, if a solution to 
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those conditions is found, the trapped surface area gives (the lower bound to) the entropy 
production in the collisions. Technically construction of the trapped surface closely follows 
early works in flat space background[15, 16, 17, 18]. The specific problem addressed in that 
work [9] was central collision of two point black holes. Among its important conclusions 
was e.g. a prediction of the entropy dependence on the (CM) collision energy S(E) ~ E 2 / 3 . 

This approach has been then generalized to the non-central collisions. We found [11] 
that trapped surface formation is not possible beyond certain critical impact parameter, 
depending on the collision energy. Furthermore, the disappearance of the trapped surface 
happens suddenly, as a 1-st order transition. An intriguing observation, also pointed out in 
our paper [11], is that phenomenologically the multiplicity of the produced particles (the 
entropy) per participant nucleon in ultrarelativistic collisions at RHIC also changes rapidly 
between "non-thermal" peripheral and "thermal" more central collisions. The specific re- 
sults about the trapped surface were later confirmed in [10, 12]. Like it has been the case in 
flat space, the value of the critical impact parameter can be understood as a bound on the 
angular momentum for the shock wave pair at a given center of mass to form a AdS-Kerr 
black hole. 

All the above-mentioned works have been using a shock wave arising from a point 
source in the bulk (small black holes) . The size of the colliding nuclei were thus incorporated 
via the distance of those objects from the boundary along the holographic coordinate z. 
However, as emphasized in [11], this is an oversimplification of the problem. The transverse 
extension of the colliding objects in x 2 ,x 3 can be introduced independently of the profile 
in the holographic z direction. The latter, due to very basic features of the AdS/CFT 
correspondence, should be ascribed instead to the intrinsic scale variable, in the sense of 
the renormalization group, describing its microscopic structure. In the collision of ordinary 
objects it would be interatomic scale, for high energy QCD the holographic coordinate z at 
which the colliding object are before the collisions should represent the typical scale of their 
wave function, known as the "saturation scale". This scale affects the typical "equilibration 
time" and other properties of the problem, thus should be taken into account. Other scales 
can also be important for the equilibration process, e.g. the longitudinal width of the 
nucleus was studied in a more recently paper [13], which contains a numerical evolution 
of the nonlinear Einstein equation and found interesting behavior of slowing down of the 
nucleus after collision. 

In this paper, we will focus on the effect of the saturation scale and model it with 
the simplest possible geometry, proposed for this purpose in [11]. It is a collision of wall 
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shock waves, which are infinite and homogeneous in 2 transverse spatial dimensions. The 
extension in z of the trapped region has been found for the collision of such wall shock waves. 
It has been done for the simplest case of a symmetric collision, in which both colliding walls 
are the same. We start in this paper discussing a more general case, in which two colliding 
walls are not the same. Physically, one may think of two colliding objects made of different 
materials with different densities, which are modeled by their different "saturation scales" 
jj- , ■ The question we will answer is the precise critical condition on their values z\ , Z2 
beyond which the trapped surface is not formed. 

Perhaps the reader may wander why are we interested in such a question. It is clear 
that one of the most important variable is the energy (rapidity) of the colliding objects: the 
black holes can only be formed if it is large enough. However, let us also remind the reader 
that in heavy ion collisions the energy per nucleon is not the only important variable: for 
example rapid equilibration and hydro dynamical behavior experimentally observed at RHIC 
for collisions of two heavy ions such as AuAu, are indeed not observed say for deuteron-Au 
collisions at the same rapidity of the colliding nuclei. Similarly, we find that two walls, 
made of sufficiently different materials, can also collide without classical equilibration and 
entropy formation, at the same energy at which the symmetric walls would produce the 
trapped surface. 

Another issue, to be addressed in section 4, deals with the difficult problem of find- 
ing the gravitational solution for the non-zero time, in the future quadrant of the time- 
longitudinal coordinates. In flat Minkowski space-time this is a long-standing problem of 
the general relativity. Recent numerical studies [19] have managed to reach gamma factor of 
the order of few units and reasonable agreement is observed with previous partial analytic 
results reported in [14, 16, 17, 18]. However, the problem gets even more complicated in the 
curved 5-dimensional space AdS§ needed for current applications, see e.g. [21, 22, 25, 13]. 
The "Next-to-Leading Order" (NLO) effect we will discuss are the "debris" produced in 
the shock wave collision, the gravitons radiated perturbatively. We will compute the NLO 
correction to the metric, and read the corresponding stress tensor on the dual field theory 
on the boundary: such "early time" stress tensor plays an important role in the theory 
of heavy ion collisions, as it provides the initial conditions for the standard hydrodynam- 
ical treatment. We will follow most closely the work by Taliotis[25] in the settings: the 
main difference is that our source is localized in the holographic direction z, instead of the 
transverse directions as in [25]. 
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2 Wall on wall shock wave collision 

We start with the wall shock wave model proposed in [11]. The metric of a single shock 
wave moving in direction x + is given by: 



, 9 r9 -dx + dx + dx 2 : + d>(z,x + )dx +2 + dz 2 . . 

ds A = L (1) 



The shock wave profile cf>(z) satisfies the following equation: 



(d 2 - *d z )<f>(z,x + ) = -WttG 5 ^5(x + )S(z - z ) (2) 

with RHS being the source, which has infinite extension in the directions of x±, thus 
the name wall shock wave. The parameter zo is interpreted as the inverse saturation scale. 
The solution to (2) is given by: 



{z,x + )=A^G b ^5{x + ){ 



2 4 

z < z 



(3) 

1 Z > Zq 



The stress tensor follows from (3) reads: 



'-++ 



»6(x + ) (4) 



Now consider the collision of two shock waves, as a model of heavy ion collisions. The 
metric of the shock waves before collision is given by: 

2 2~ dx + dx~ +dx 2 , + 4>\(x + , z)dx +2 + (f)2{x~ , z)dx~ 2 + dz 2 
ds = L k (5) 

z z 

The shock wave profiles solve the following equations: 



(d 2 - h z )M* + , z) = -167rG 5W ^<5(x + ) ( 5(z - Zl ) (6) 
(8 2 - ^d z )<h{x + , z) = -l6irG 5 ti 2 jts(x-)5(z - z 2 ) (7) 

Note we have absorbed the delta function into the definition of the shock wave profiles 
as in [25]. The dual stress tensor reads: 
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T__ = /i 2 <5(x-) 



(8) 



The superposition of two shock waves (5), solves the Einstein equation in the region 
with 9{x + )9{x~) = 0. (Here 9{x) is the Heaviside step function, 1 for positive and for 
negative argument.) The shock waves only interact and modifies the metric in the future 
quadrant 9{x + )9{x~) > 0. 

As explained in the Introduction, our colliding walls are dual to "nuclei" of infinite 
size, so the concept of impact parameter does not exist. Instead, we have also chosen two 
nucleus to have the same energy fj,± = (12, but different saturation scales Z\ 7^ z-i- (To be 
specific, we demand z\ > Z2-) Although our shock waves are sourced by the delta functions, 
they have finite size in the z direction, decreasing both into the UV and the IR. This is 
different from other approaches using sourceless shock waves[20, 21, 22, 23, 24, 13]. 

These finite extension of the shock waves in z explains why the trapped surface can 
be found also in a finite interval in z, we will call upper and lower positions of the trapped 
surface z a , Zb- 

The entropy lower bound, dual to the "area" of the trapped surface is given by: 



3 Critical condition for trapped surface formation 

In this section, we will construct the trapped surface associated with the collision of two 
shock waves. Let us for the completeness recall the mathematical basis defining the trapped 
surface. The equations are produced by required vanishing of the so called "expansion" com- 
bination: loosely speaking it means that the geodesies of forward moving, outgoing massless 
particles should converge on this surface. The limiting case when the geodesies neither con- 
verge nor diverge defines the maginally trapped surface. It can be shown to correspond 
to a relatively simple problem a la electrostatic solution in a cavity (the Laplacian with 
given sources) with zero boundary condition on the surface, complemented by additional 
nontrivial condition for the magnitude of the field derivatives at the surface itself. Following 
[9, 11, 10], the master equation for trapped surface is given by: 



S = 



2A f ^/gdzd 2 x± 



4G 5 2G 5 



s = 




(9) 
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zH'l - z% - 3*, = -\^G hiH z\b{z - Zi ) 

^(z a )%(z a )^ = ^(^(^l = 4 (10) 

with i = 1, 2. The trapped surface for wall-on-wall shock wave collision is just z a < z < z b . 
The first two equations can be solved as: 



Ci ( 4 - ) Z < Z,; 



with 

" 1)^6 

^ = -47rC 5 ^ " 4 _ (12) 

z b fa 
(ff " l)*o 

A = -47rC 5 ^ z ° (13) 

fi fa 

We can always apply a longitudinal boost such that both shock waves have the same 
energy density y/niH2- Then the third equation in (10) leads to 



dC 2 = D 1 D 2 = ^ (14) 



Let us consider the case z\ = z 2 = zq first. (14) leads to: 



z a + z h = = A 1 (15) 

(z a + z b ) 2 - 3z a z b _ L 3 _ 1 

(^ 6 ) 3 4 A 2 1 ° j 

in which two appearing combinations of parameters are for brevity called A±,A 2 . The 
resulting cubic eqn 

to) 3 + 3A 2 to) - A X A 2 = (17) 

can be solved by Cardano formula. The explicit solution is not illustrative and is not showed 
here. We note, however the solution has to satisfy the inequality 4:Z a z < (z a + z ) 2 = A\, 
which gives rise to the following constraint: 
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2ir 2 , 

-J^l^l > 1 (18) 



where we have used G5 = jj^. (18) is the critical condition for trapped surface formation 
in a symmetric collision of gravitational shock waves. 

z 4 z 4 

When z\ > z 2 , we define = M, -^1 = A2, (14) can be simplified to: 



z a z b Z a z b 
Za\ 2 I f 2 _(_ 1 _ A1+A2 + I 

z b "•" \z a A1A2 



L 2 (19) 



(87rG 5 M) 2 (l-AiA 2 ) 



where the first equation follows from = 1 and the second equation can be obtained 

from dC 2 = The first equation can be used to give f^ + f^ = ^ (Al+ ^ 2+ ^. Combining 
this with the second equation, we can express /iasa function of Ai and A2, which in terms 
of variable F = AiA 2 and r = {X ^lf = ^feF reads: 



( wt^)» _ w_ ^1 + ,j 1/2 ^ _ 3 ) (2 „, 

Note that r depends on the degree of the asymmetry of the collision, for z\ = z 2 one has 
r = 4. Let us thus fix r and study the RHS of the (20) as a function of the other variable F, 
to be called A(r,F). Prom the second equation of (19), we know F < 1 and by definition 

3 

F > 0. For a given r, we have the following limits: as F — > 0, A — > F~i and as F — > 1, 
A — )■ Yr7?(v / r : + 2) 2(\/r — 2). Unless r = 4 (the symmetric case), in both limits the function 
tends to positive infinity. Therefore a minimum must exist at certain F = F m i n , which gives 
rise to the critical condition we are looking for. Fig.l contains a plot of A as a function of 
F at several r. 

The extremum of A(r, F) is found to be the roots of the following equation: 

-3F - 4rF 2 + 3VrF + 3F 3 + 13F 2 - GFVrF + 3F 2 \/^F + 3 = (21) 

It is not difficult to locate the minimum of A(r, F) numerically, which gives rise to a critical 
condition for the collision energy: 

^M^) 3/2 > VW) (22) 

where we have used G5 = G(r) is the minimum of A(r,F) at a given r. For 

r = 4(z± = z 2 ), A(r,F) has a minimum at F — > 1: G(r) = 4. We recover the critical 
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5000 




r=4 ■ 



r=6 



1=10 ■ 



i=20 



Figure 1: A as a function of F. A minimum always exists in < F < 1 for r > 4. In the 
extreme case r = 4, the minimum locates at F = 1 



condition for the symmetric collision (18). For general r > 4, we find G(r) numerically and 
as grows, \J G{r) has a power like asymptotics \J G(r) ~ f f 1 ) • Fig. 2 shows a the 



power law dependence of G(r) on ^. The power measured by the slope in the log-log plot 
is approximately 1.5. 

The asymptotic power law behavior of G(r) can be obtained analytically. We note 
the root of (21) corresponding to the minimum of A(r,F) goes to zero as r — > oo. As the 
result, (21) simplifies to — 4rF 2 + 3\/rF = 0, which is solved by 



3 

ri + 



F 



3\ 2/3 i 



(23) 



where ••• denotes subleading terms. Substituting the root to A(r,F), we obtain G{r) 
. Combined with the definition of r, we indeed have: 



y/Gfr 



3 



z 2, 



In the limit r S> l(z± S> Z2), the critical condition simplifies to 

16tt 2 



•n4 > 1 



(24) 



(25) 



We would like to point out the non-uniqueness of the trapped surface, as first re- 
marked by Eardley and Giddings [16], the unusual boundary value problem defining the 
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Figure 2: A log- log plot of y/G{r) versus | 



trapped surface could have multiple solutions. We will see it is indeed the case in our 
wall-on- wall collision.* 1 . 

Suppose we have the energy of the shock wave well above the critical valuea.e. Aq = 



^Tli-{z\Z2^l 2 » G(r). We know from the previous analysis that 



A(r,F) -»• F" 3/4 as F^O 

A(r,F) — >■ -— !— -(a/t + 2)5 (-y/F - 2) as F ->1 

1 — X* 

This allows two solutions F = A~ 4/3 + • • • and F = 1 - (y^My^zg) + _ . . _ Wit hout 
explicit solution of the trapped surface, we can compare the area of two corresponding 
trapped surface, which is related to the entropy production per transverse areafll]: 



27T 1*2 Z} 



b ■ 

1/2 



= (27riV c V) 1/3 A~ 1/3 F- 1 / 4 (l + v^F - 3F) (26) 

With the former solution, we have s = {Q.TrN^fij 1 ^ 3 + • • • , while the latter solution 
gives rise to s = (2vriV c V) 1/3 A) 1/3 (\ /F - 2 ) 1/2 + • • • • In the limit -> oo, the former 
trapped surface has a much greater area than the latter. Therefore we choose the former as 



#1 Apart from this, there is also the foliation dependence of the trapped surface, which we do not discuss 
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Figure 3: (color online)A view of the outer-most trapped surface formation in a wall-on-wall 
collision with two sources at different depths. The pink and blue area indicate the growth 
of the trapped surface \E' in the bulk. The sources of the shock waves lie at ^ = 1.3 and 
^ = 1.6, and the energy density is fixed by ( 87rG 5A^( z i^) _ 20. The trapped surface at 
the collision point is bounded by z a and Zf, 

the "outermost" trapped surface. This branch of solution is precisely the one used in [26] 
for a comparison of sourced shock wave and source-free shock wave. 

A view of the trapped surface formation is included in Fig. 3. We note the trapped 
surface only starts to appear when the wave fronts of the shocks are separated by finite 
distance, unlike the situation in point shock wave collision, where the trapped surface starts 
to grow even for infinitely separated shock waves. Furthermore, the origin of the critical 
condition in wall-on-wall collision is different from that of point shock wave collision : The 
latter can be understood as a constraint on angular momentum for the shock wave pair 
to form a AdS-Kerr black hole. The energy dependence of the critical impact parameter 
obeys an asymptotic power law, with the power extracted numerically in [11] to be 0.37, 
and argued by Gubser, Pufu and Yarom to be 1/3[10]. The power 1/3 was confirmed later 
in a detailed numerical analysis [12]. The critical condition for wall-on- wall collision (22) 
has no analogy here, as the shock wave pair does not have an obvious angular momentum. 
In the next section, we will compute the NLO stress tensor in the dual field theory, which 
will help us to understand the origin of the critical condition. 
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4 The NLO stress tensor after the collision 



The Einstein equation in the presence of the cosmological constant is given by: 

GW + 6<v = -8vrG 5 «V (27) 

with the 5-dimensional source in the bulk. We choose to work with an alternative 
form of the Einstein equation: 

iV - 4 9»v = -8ttG 5 (28) 

where S^ v = J^ v — \ Jg^ u - We have set the AdS radius L = 1. Since the relevant scales 
are completely fixed by \x and Zi, we expect L will be absent in the final result of the stress 
tensor. The equation (28) to the first order in the amplitude of the shock wave is: 

R$ ~ 4$ = -SttG^J = -8vrG 5 (J« - i jWg$) (29) 

where the upper index denotes the order of the quantity with respect to the amplitude of 
the shock wave. e.g. gffi is the pure AdS metric. Before the collision the superposition of 
two shock waves are solution to this equation. Their contribution to the source is of first 
order: 

&kG 5 J$. = -\v 2 4>i 

8itG 5 J { }1 = -^V 2 02 (30) 

After the shock waves pass through each other, the source of either shock wave feels 
the field of the other shock wave and deviates from its original trajectory. The deviation 

(2) 

gives rise to the second order correction to the source: J^J. On the other hand, since the 
superposition of two shock waves does not satisfy Einstein equation after the collision, thus 
a nonvanishing Ricci tensor is expected from the superposition. We collect this contribution 
into R^i) 1 ^, which can be interpreted as the interaction of the shock wave fields. Combing 
the two contributions, the second order Einstein equation now takes the following form: 

+ - 4g$ = -8nG 5 S$ = -8*G b {J$ - _ jWgW) (3 1) 
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Our contracted source of order k is always denned as = jj$ g^i ) . 

(2) 

The calculation of J^J needs some explanations. Since our wall shock wave has trivial 
dependence on transverse coordinates, the problem gets simplified a lot. The sources can 
only move in z, and the determination of their trajectory after the collision is subtle. In 
general it depends on the equation of state of the extended source itself. We will assume 
the action due to the shock wave source is of Nambu-Goto type, which is proportional to 
the invariant of the extended source. 

Calculating the geodesic of one shock wave source in the background of the other 
shock wave /zeZeZ(details can be found in the appendix). The second order source is given 
by: 

8ttG 5 J { +1 = ^(j fodx-d+V 2 ^ + ^jdx-j dx-dzfodzV 2 ^) 
8ttG 5 J {2) _ = ±(J ^dx+d-V 2 ^ + \j ' dx+ J dx+ d^d.V 2 <f> 2 ) 
8irG 5 j| 2 2 = ^(0iV 2 </> 2 + 2 V 2 0i) 



8^G 5 j| 2 j = ±J dzfodx-V 2 ^ 
r( 2 ) _ 1 f a a.a„-ti1. 



8ttG 5 JT£ = -J d z <f)idx-V 2 fo (32) 
We can check the following relations to the second order: 



(v^j„)( 2 ) = v^)j( 2 )+v^)j« 

= -2z 4 5(z - Zl )S(z - z 2 )(8(x + )5(x-) + 8( X -)5(x + )) = (33) 

or.w (2) = <r {0) j$ + <r (1) = o (34) 

The first one is the conservation of the source, which is a necessary condition for the 
consistency of Einstein equation. The second traceless condition allows us to simplify the 
RHS of (31). Moving The Ricci tensor quadratic in the first order metric to the RHS 

and noting the tracelessness of = = 0, we obtain the reshuffled Einstein equation 
for the second order corrections only 

R$ ~ ±9$ = -&xG 5 J$ - = -&kG 5 J$ (35) 

where we have defined the effective source jffi = jf$ + 8 J~ G Rfy 1 ^ 

It is easy to work out R^) 1 ^ for the case of wall shock waves and we obtain the effective 
source as: 
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8vrG 5 J| 2 | = dx-fod+V 2 ^ - i J f dx~ J dx~ d z <fod z V 2 </>i 
8ttG 5 J {2) _ = ^J dx+0i5_V 2 2 - ^ J dx+ J dx+dzfadzV 2 ^ 

8ttG 5 j^_ = h<h.v 2 <h + <A2V 2 0i) - d+fad-fa + d z fad z fo - -(<hd z <h + <h%fa) 



SvrGsJg = \ I dx-d z <t> 2 V 2 <t>i ()-<>i<)-o-2 



8vrG 5 J {2) z = \J dx+d^V 2 ^ - d-<t> 2 d z <t>± 
8vrG 5 J| 2 | = ^(0i9 2 2 + <hd x <h) 

8vrG 5 Jf) = -2(0 2 a>i + fc&M + *{<hd z <h + <hdz<h) ~ 2d z <hd z <h (36) 

From here on, we can use the method developed in [27, 28] to compute the stress tensor 
on the boundary field theory to the NLO. The procedure is to first obtain the reshuffled 
source s mn defined as 

s rnn = ^mn ~ J {^zm,n Jzn,rnJ + '^ l ,m,n + ^Imti^z (37) 

with 

h = \J*dz-z(j®- v mn J£l + * £ (-V mn Jz 2 2,n) ) (38) 
In our particular case, h is given by: 

8nG 5 h = -401^2 + 4 J dz-z J dzhd z fad g <j>2 (39) 
The reshuffled source takes the following form: 

8-kG^,s ++ = i j dx~ 4>2d+V 2 4>\ + ^ J dx~ j dx~ d z 4>2d z V 2 4>i 

— J dx~ j dzd z (f>2d + V 2 (f>i — 2 J dz&dzd+fa + 2 j dz-z j dz^-d\d z 4>id z 4>2 
8vrG 5 s__ = i j dx + (pid-V 2 <j)2 + dx + j dx + d z (f>id z V 2 (f>2 

- J dx + j dzd z (f> 1 d-V 2 (f>2-2 J dz(f)id z d 2 _<t)2 + 2 j dz ■ z J dz^d^fod^i 
87rG 5 s+- = difadzfa - 2d+cf> 1 d-<f>2 + \ J dzfad^fa + <hd g v 2 <h.) 

+2 J dz ■ z J dz^-d + d z <pid^d z (/)2 - J dz^-d z (f>id z <p2 

8irG 5 s ±± = 2 J dz-d z ^ x d z ^2 (40) 
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Separating the derivatives of x + , x and derivative of z, we obtain: 
87rG 5 s ++ = ^4>2^ 2 4>i0{x')5'{x + ) + ^V 2 ^- 6(x-)5(x + ) -2 J dz&faS" \x + )S( X -) 
+2 j dz-z j dz- z ^ 2 5"{x + )5{x-)- j dz4> 2 V 2 fa9(x-)5'(x + ) 

8ttG 5 S— = ]^)iV 2 foO{x + )5'{x-) + ^' l V 2 $ 2 x + 9{x + )5{x~) -2 j dz$' 2 fa6" \x~)5(x + ) 
+2 f dz-z f dz-4>' 1 4>' 2 5"(x-)5(x + )- [ dz$\V 2 faO{x + )&{xT) 



8vrG 5 s + - = ^ dz^V 2 ^ + 4> 2 V 2 ^ 1 )5(x + )5(x-) + ^' 2 5(x + )5(x-) 
-20 1 2 S'(x + )5'(x~) + 2 Jdz-z J ^i^(x + )S'(x-)- J dz- z ^' 2 5{x + )5{x 



1 



8ttG 5 s ±± = 2 / dz-4>' 1 4> 2 5(x + )5(x-) 



(41) 



where (pi(x + ,z) = (f>i(z)5{x + ) and 4>2{x ,z) = 4>2(z)8(x ). In (41), all primes are 
ordinary derivatives. The explicit forms of </>i and 4>2 are given by: 



4>i = 4irG 5 m 



b 2 = 47rG 5 /i 2 



z\ z 4 - (z 4 - z 4 )9(z - zi) 

I? zf 

4 z 4 - {z 4 - z 4 )9(z - z 2 ) 



(42) 
(43) 



The reshuffled source (41) will be convoluted with a bulk to boundary propagator in 
AdS. Such propagator has been built in various applications of AdS/CFT, e.g. [29, 22, 28, 
30, 32]. Propagators in an AdS shock wave background were found in [33, 34]. We will use 
a slightly different propagator from the above. The propagator takes the following form in 
the lightcone coordinates: 



Pr 



9(x^ 



x +/ + x~ 



x-') 



2tt 



8"'{z 



8w 3 



w) 1 / 2 ) 36" {z 



w) 



8w 4 



+ 



W(z 



w 



8w 5 



(44) 



with w = \/(x + — x + '){x~ — x~') — (xj_ — x' ± ) 2 . 

The details of the propagator are included in the appendix. Since we are dealing with 
wall sources, which do not depend on x' ± , we can perform the integral with respect to the 
transverse coordinate x' ± . By repeated use of integration by parts, we end up with a concise 
form: 
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//"CO 
d 2 x' ± P R = 2vr J dx ± ■ x' ± P R = 6{x + - x + ' + x~ - x~') x 

S'(Z- y/(x+ -X+')(X- -X-')) S"(Z- yj{x+ -X+'){X- - X~')) 



(45) 



8(x+ - x+')(x- - x-'f' 2 8(x+ - x+')(x- - x-') 

The final task is to convolute the source (41) with the integrated propagator (45). 
Due to the presence of the delta function, the integration in z is trivial. We are only left 
with integration of x +l and x~'. Completing the integrals, we obtain as the final results: 



pNLO = 8TT 2 G 5 flifl 2 
pNLO = 8TT 2 G 5 flifl 2 

TV? 



- x -^(z 2 -t)9(t) + ^5(z 2 -t 

-x +2 e(z 2 - t)9{t) + X+Z2 



2x~ 



S(z 2 



NLO = }_ T NLO _ 87r2 C5^1^2 



2t 2 0(z 2 - t)9(t) - ^5{z 2 - t) 



(46) 



+- 2 -±± N 2 

where r = Vx + x~ is the proper time. This is the main result of this chapter. 

(Few technical comments on the derivation: We have also used the distributional 
relations in the final results S^ n \x)f(x) = (—l) n f( n \x)5(x). In doing this, we have treated 
r and 2p as separate variables. It is however necessary to keep in mind one subtlety. We 
have assumed the source has a series expansion near the boundary z = 0, in the derivation 
of the propagator. As our source contains delta functions and Heaviside theta functions, 
(46) is obtained with a particular representation of them and the limit is taken in the final 
results.) 

Several more general comments on the result, the NLO stress tensor (46), are in order: 

i) The NLO stress tensor is conserved and traceless d m T^° = 0, r\ mn T^° = 0. The 
presence of the delta function is necessary for the conservation relation. As conjectured in 
[26], The limit z 2 — > oo of our results should recover the NLO stress tensor in the collision of 
sourceless shock wave [21, 22]. We can see it is indeed the case as 0{z 2 — t) = 1, S(z 2 — r) = 0. 

ii) (46) is actually boost invariant. In a comoving frame with coordinate r = Vx + x~ 
and r\ = ^ In ^p, the NLO stress tensor takes the following form: 



T^o = 8 -^l^ [ 2r 2 e(z 2 -r)e(r)] 



T „ LO = SfG^ _ t)6{t) + 2t24s{z2 _ r)] 



Tff° = *^ G >gV» [Ar 2 6(z 2 - r)6(r) - z\^ 2 - r) 

^ c 



(47) 
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The boost invariance is a special property of the NLO stress tensor, which is symmet- 
ric under the exchange of the two shock waves. Since we are colliding asymmetric nucleus, 
we expect higher order correction should violate boost invariance. 

hi) It is interesting to note that the NLO stress tensor does not depend on z\. It 
suggests the NLO stress tensor for collision of two nucleus with different saturation scales 
does not feel the softer saturation scale 

iv) The appearance of the Heaviside theta function is of particular interest. It encodes 
information on thermalization. As the LO stress tensor = Hi8(x + ), T LO = n 2 6{x-) 
has a simple interpretation as nucleus moving on the lightcone. The NLO stress tensor (46) 
tells us matter created in the collision is only nonvanishing when < r < z 2 . At time t > z 2 , 
matter created in the collision separates into two pieces z 2 < X3 < t and — t < X3 < — z 2 . 
Presumably higher order correction is needed to fill the gap. This also suggest the NLO 
result is insufficient to provide an initial condition for hydrodynamics. 

v) Comparing the normalization of the delta functions in the LO and NLO stress 
tensor, we conclude the perturbation should break down when fi < SttG^/j^z^, which is 
precisely the critical condition (25). 

Therefore the field theory interpretation of the thermalization condition is understood 
as the breaking down of perturbative treatment. Presumably the combined effect of all the 
gravitons should be included in further evolution of the trapped surface, from its position 
at time zero discussed at the beginning of the paper. 

Alternatively, we can take a bulk point of view: The perturbation breaks down when 
the sources of the shock wave, originally moving at constant radial position, deviate signif- 
icantly in the radial direction. Specially, we have worked this out in the Appendix A. With 
A*i = A*2 = A 1 ) the sources of shock wave gain velocities 



after the collision. Due to the special profile of the shock waves, the source deeper in the 
bulk does not shift its path in the NLO computation. The perturbation breaks down when 
u 2 ~ 1> which again is consistent with the critical condition (25) when |^ S> 1. 

The improved understanding of the critical condition leads to the following prediction: 
In the collision of two nucleus with the same energy density but different saturation scale, 
the thermalization condition is insensitive to the softer saturation scale. The energy density 




(48) 
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has to exceed certain critical value set by the harder saturation scale as (25) in order to 
reach thermalization. 

5 Discussion 

In this work, we have constructed the trapped surface in a wall-on-wall collision, which 
is used to model collisions of nucleus with different saturation scales. We have derived a 
critical condition for matter equilibration in nucleus collisions. The condition (22) is set 
by the saturation scales of both nucleus. The critical energy scales as the ratio of the 
saturation scales approximately by a power law, with the power 3/2. The approximate 
power law indicates the critical energy is insensitive to the softer saturation scale. We have 
also observed a non-uniqueness of the trapped surface when the energy density is beyond 
the critical value. The outer-most trapped surface is selected for an estimate of the entropy 
production. 

We have computed the NLO stress tensor on the boundary. The result turns out to 
be independent on the soft saturation scale l/z\. Based on the NLO results, we propose the 
critical condition corresponds to the breaking down of the perturbation, i.e. when the LO 
and NLO correction become comparable. The criterion reproduces the critical condition 
(25). On the other hand, the critical condition is also understood in terms of bulk physics. 
The breaking down of the perturbation is encoded in the condition when the sources of the 
shock wave gain significant deviation in its velocity after the collision. This also leads to 
the correct critical condition (25). While in the NLO, no dependence on z\ is observed, 
it must show up beyond NLO, as the source deeper in the bulk will also deviate from its 
original path, giving rise to correction to (25). It is tempting to see how this shows up in 
higher order computation. 

Finally we stress the physics of critical condition for matter equilibration is very 
different from the counterpart in the collision of point source shock wave. In terms of the 
gravity dual, the latter originates from the constraint on the angular momentum possessed 
by a pair of black holes in order for the merging to be possible. The critical condition 
for wall shock wave collision can be understood as the breaking down of the perturbative 
calculation. In the dual field theory, it manifests as a constraint on the collision energy for 
given parton saturation scale. 
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A NLO source in shock wave collision 
A.l Point source shock wave in AdS 3 

It is helpful to look at collision of point shock wave in AdS% first. The LO metric is given 
by: 

2 _ -dx + dx~ + dz 2 + (f>idx +2 + fodx' 2 , 
ds — 7; (49 J 

z z 

The shock wave profiles <p\ and <p2 are normalized as: 

V 2 0i = -16ttG 5 S(x + )S{z - zi) (50) 
V 2 2 = -16ttG 5 5{x-)S(z - z 2 ) (51) 

where V 2 = d 2 — \d z is a Laplacian operator. The NLO source arises from the 
deviation of the path of one shock wave source in the presence of the other. For point 
source, the null geodesic equation is given by: 

^ + Y^vP = (52) 

For source of shock wave 1 before the collision, x~ can be chosen as the affine param- 
eter A, thus u~ = 1, u + = u z = 0. After the collision, the geodesic to the first order in the 
shock wave amplitude is as follows: 



Assuming A = holds after the collision, we find from (53) and (54) that u + = <p2 
and u z = f dx~^p-. However we see it contradicts (55) as = 0. This indicates that 
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x~ is no longer a good affine parameter, but to the order we are interested, u + = 4>2 and 
u z = f ^p^dx~ remains valid, as correction will be of higher order. Integrating once, we 
further obtain x + = f dx~4>2 and z = Z\ + / dx~ f dx'^j 1 . 

The covariant source due to shock wave 1 has the general form: 



juv = # u » u »S(x + - X + {x-))5{z - Z(x-)) (56) 

where X + (x~) and Z(x~) specifies the trajectory of the point source. # can be some 
function of x + , x~ and z. Writing the LO covariant stress tensor is simply: 



SnG 5 J—^ = -2z 4 V 2 0i (57) 
8itG 5 J ++{1) = -2z 2 V 2 (p2 (58) 

The NLO source comes from the correction to and x^. Adding the contributions 
from two shock waves, we obtain: 

8vrG 5 J— (2) = 2z 4 ( J dx-fcd+V 2 ^! + dx~ J dx~ d z fod z V 2 fa) 

8ttG 5 J ++(2) = 2z A { J dx + fad-V 2 fa + ^Jdx+J dx + d z fad z V 2 fa) 
8irG 5 J + -W = -2z\faV 2 fa + faV 2 fa) 
8itG 5 J- z W = -2z 4 ■ 1 J d X -d z faV 2 fa 

8ttG 5 J +z(2) = -2z 4 • 1 J dx+d z faV 2 fa (59) 
The conservation of the source to the second order can be checked (V M T^)( 2 ) = 

q ^m^(2) _|_ y^t Xu ^ + p^t^ 2 ) + p^r^ 1 ) + y u ^t^ x ^ = o 

M fi\ fj,X (i\ fi\ 

A. 2 Wall source shock wave in AdS$ 

Now we look at wall shock wave in AdS$. The LO metric is given by: 

, 2 — dx + dx~ + dx 2 , + dz 2 + fadx +2 + fadx~ 2 . , 

ds = k (60) 

z z 

The shock wave profiles fa and fa are normalized as 

V 2 fa = -16ttG 5 6(x + )S{z - zi) (61) 
V 2 2 = -16ttG 5 5{x-)5(z - z 2 ) (62) 
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The Laplacian operator becomes V = d\ — \d z due to the additional transverse 
directions. Being different from the point source, the trajectory of the source is specified 
by X fl (a) with a the worldvolume parameters. The induced metric is given by: 



af3 ~ da* da^ v 



I —u + u +u z 

V 



(63) 



with v? 1 = ^jj- and u z = 4|. Assuming the action of the shock wave depends on 
det h only, then the trajectory can be effectively determined by considering a point source 
in the metric 



ds 2 



-dx + dx + dz 2 + 4>idx +2 + <p 2 dx 2 



(64) 



Working out the geodesic deviation, we find surprisingly that the trajectory of the 
wall source is the same as point source in AdS%. As a result, the LO and NLO source are 
given by: 



8^G 5 J— (1) = -2z 4 V 2 ^i (65) 
8^G 5 J ++(1) = -2z 2 V 2 4> 2 (66) 
8ttG 5 J— (2) = 2z4 (J dx-fcd+V 2 ^! + dx~ J dx~ d z fcd z V 2 <pi) (67) 

8^G 5 J ++(2) = 2z 4 ( J dx+fad-V 2 ^ + ^Jdx + J dx+d^d.V 2 ^) (68) 
8^G 5 J + ~ (2) = -2z 4 (^V 2 2 + 2 V 2 0i) (69) 
8ttG 5 J- z{2) = -2z A ■ ^ J dx~d z <j) 2 V 2 (f>i (70) 

8^G 5 J +z{2) = -2z 4 • 1 J dx+dzfaV 2 ^ (71) 

While (65) has the same functional form as (59), they are different in the Laplacian 
operator. We can check (65) is again conserved and the Christoffels involving the additional 
directions are accounted for the difference in the Laplacian operators. 

With some care, we can obtain the NLO contravariant source, which include contri- 
bution from both LO and NLO covariant sources. The result is shown in (32) in the main 
text. 
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B The bulk to boundary propagator 



In this appendix, we want to build a propagator, which produces the stress tensor on the 
boundary field theory when convoluted with the bulk source. We start with a bulk to bulk 
propagator for massive scalar defined as follows: 



^dJV^gg^d^G - m 2 G = -^5 {d+1) (x - x') 



(72) 



The metric is the Poincare patch of AdSd+i ■ Using the Fourier transform: G(u, k, z) 
f G{t,x,z)e- iul( - t - t ^ +i ^ £ -^\ (72) takes the following explicit form: 



z 2 (u 2 - k 2 )G + z 2 d 2 G + (1 - d)zd z G = m 2 G = z' d+1 5{z - z') (73) 

The boundary condition to impose is that G — > as z — > and G is outgoing as 
z — > oo. The solution to (73) is found to be 



G = -{zz')h A (Vk 2 - u 2 Z< )K A {^k 2 - uj 2 Z> ) (74) 

where A = ^ d2 + 4m2 an d z> = max{z, z'}, z< = min{z, z'}. 

The inverse Fourier transform gives the bulk to bulk propagator: 



/ \ — 

G(t,x,z) = -feb- I lA(Vk 2 - ^z^K^k 2 - u; 2 z > )e tuJ ^- i ^-" ; 'Uu;d d - 1 k (75) 
(2vr) a J 

Note there are two branch cuts on the real axis (— oo, —k) and (k,oo). The retarded 
propagator can be obtained if we take the integration contour of u slightly below the real 
axis: ui — > ui — ie. We can push the integration contour to wrap around the two branch 
cuts, so that all the contributions come from two sides of the branch cuts. 



G(t,x,z,z?) = -j^fO(t-t')( [ dui+ dw) [ d d-i ke i"(t-t')-ik(3-S') x 

l^ 71 ") J -oo Jk J 

K A (-iVk 2 - w 2 z>)/ A (-i\/w 2 - k 2 Z< ) - K A (iVk 2 - uj 2 Z> )I A (iVoJ 2 - k 2 Z< )) 

= - i 4f^9(t-t') r du [ d d ~ 1 k2iTj A (Vuj 2 - k 2 Z> )J A {^/u 2 - k 2 Z< ) x 
(2vr) a J k J 

sinw(i - t>) e -ik2-*') (76) 
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Doing the angular integration for the spatial momentum k, we obtain: 



G(t, x, z, z') = - 0(t - t') dco [ k d - 2 dkdn d - 2 sin oj(t - t')e ikrcose x 

(2vr) a! 1 J k J 

J A {Vu 2 - k 2 Z> )J A {VoJ 2 - k 2 Z< ) 

= - -^0^ O(t-t') J™ du J k d ~ 2 dkde(sme) d ~ 3 dn d - 2 smu;(t-t')e ikrcose x 

Ja{V^ 2 - k 2 Z> )J A {Vu 2 - k 2 Z< ) 

(zz')i 2^71^ f°° f d-i 

= - , J_! — A — 9(t-t'-r) / doj dkk~ J d-3_(kr) sin u(t - t') x 

(27TJ f~2~ Jk J 2 

Ja{Voj 2 - k 2 z> )J A {Vu 2 - k 2 z< ) (77) 
We have defined r = \x — x'\. Writing /3 = \J uj 2 — k 2 allows us to do the /c-integral[35]: 

d d-1 d-1 



G(t,x,z,z') = - ,[ ZZ 2\ 2 \V e(t-t>) [ S ^^±Et^l pd(3k d ^dk x 
J A (/3z> ) J A {(3z< ) J d-3 (kr) 



d 

( zz ' 



il-2^7r% -t'-r) j J A {pz > )J A (p Z< )plj_i_2(pw)w-^ 1 (78) 



(2vr 

where w = y/(t — t') 2 — r 2 . The final integration of /3 can also be done [35], we end 
up with 



G(t -t',x-x', z, z>) = -^L 2 ^nh(t - t> - r) x 



' \\ ZZ !) 2(sinhu) "2' sin[(— | + 1 — A)7r]e ! " 2 ,r( 3 A 2 _i(cosliu) w > Z> + Z< 

d-1 ±r! 
, . , , , —P 2 

/2^ v > \^"> A- ^ 

otherwise 



j==(zz') 2(sint;) 2 P ^(cosw) Z> — Z< < w < Z> + Z< 



(79) 

(79) is in agreement with early results on bulk to bulk propagator [31, 25]. However 
there is a non-analyticity at w = Z> + Z< , which is hidden in (79). Integration across the 
non-analyticity can lead to finite contribution, thus we choose to start with (78) in building 
the bulk to boundary propagator. 

The relevant Green's function G b (t — t',x — x', z, z') is given by: 



2 

Z — {-d 2 + d 2 x + d 2 )G b + Z -d z G h - 4G b = 5(z - z')5 d {x - x') (80) 
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The metric perturbation in the axial gauge h mn is related to the reshuffled source 
Smn by: 



(81) 



We have suppressed the tensor indices in h mn and s mn . G b is related to the bulk to 
bulk propagator by: 



G b — 2 n G\ 



|A=2,(2=4 



(82) 



Let us suppose the source adopts the following expansion near the boundary. 



s(t',x',z') = J2 s n(t',x')z' n (83) 

n 

We can perform the integrations first in z' and then in j3 to obtain: 



h(t, x, z) = - ^- J dt'd 3 x' 



6 2 n(n — 2)(n — 4) 1-n 3 - n z 2 , , 
w n b ^_K j\ LF(——,——-S;-^)8 n (lf,a/) ( 84 ) 

— ■ Z. TV 



We are interested in the coefficient of z 2 , which encodes the boundary stress tensor. 



Note lim 2 ^ ^T 1 ' 3 > ~^ 1- Tne coefficient is given by: 



-±Jd^X>Y, W n - 6 n{n - 2 l {n - 4) Sn(t<, X > 



,_ 6 n(n-2)(n-4) 1 



= -— [ dWx'Yw 

= - J- I dt'd^'dz'y w ^ n{n ~ 2) } n ~ 4) { ~ l T s{t\ x>, z')5^\z') 

Z7r / 8 n! 

n 

We can sum the n-series and obtain as our bulk to boundary propagator 



(85) 



9(t-t' -\x-x'\) 



2vr 

_e(t-t' - \ x -x'\) 

9(t - t' - \x - x'\ 
2tt 



n(n - 2)(n - 4) (-1)™ 
^ 8 n\ y 1 



J2 w n- ^ ~ l )^ n ~ 2 ) ~ 3n (n - 1) + 3n (-If 6{n) (/) 



n! 



3w" 4 



3itJ 



(5 (z — ttj) H — 8 (z — w) H — J (z — w) 



(86) 
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We can further use the property of delta function to replace 9(t — t' — \x — x'\) by 
9(t-t'): 



R _0(t-t' 



2 



-—6 [z -w) — <5 (z -w) — o [z - w) 



(87) 
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